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We estimate the isospin breaking contributions to the baryon masses which we analyzed recently 
using a loop expansion in the heavy baryon approximation to chiral effective field theory. To one 
loop, the isospin breaking corrections come from the effects of the d, u quark mass difference, the 
Coulomb and magnetic moment interactions, and effective point interactions attributable to color- 
magnetic effects. The addition of the first meson loop corrections introduces new structure. We 
estimate the resulting low-energy, long-range contributions to the mass splittings by regularizing 
the loop integrals using connections to dynamical models for finite-size baryons. We find that the 
resulting contributions to the isospin breaking corrections are of the right general size, have the 
correct sign pattern, and agree with the experimental values within the margin of error. 
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1—5 , I. INTRODUCTION 

o , 

CO . In a recent paper [3| , we analyzed the structure of the electromagnetic contributions to the mass splittings within 
isospin multiplets in the baryon octet and decuplet. In particular, we studied the leading isospin breaking (IB) 
, contributions to the mass differences that come from the Coulomb and magnetic interactions, quark mass differences, 
^ ■ and the one-loop mesonic corrections to those interactions. 

[ Our analysis was based on the standard chiral Lagrangian in the heavy baryon approximation (HBA) to chiral 

■ effective field theory as developed by Georgi Q and Jenkins and Manohar [^]. However, we used a spin- and flavor- 
Q^ [ index or "quark" representation of the effective octet and decuplet baryon fields and the electromagnetic and mesonic 

interactions rather than the usual matrix expressions for the fields. The connection of this representation to the 
usual effective-field methods was discussed in detail in 0, 0] • It has been applied in past analyses of the splittings 
, between baryon isospin multiplets [1, H, @] j and of the baryon magnetic moments [1, |7| ■ The results in each case can 

■ be summarized in terms of a set of effective interactions that have the appearance of interactions between quarks in 
Oh, the familiar semirelativistic or nonrelativistic quark models for the baryons. The results are in fact completely general 
^ in their representation of the relativistic heavy-baryon effective fields and their interactions as shown in [3, . In 

' the present case, our results connect directly to Morpurgo's general parametrization of the electromagnetic effects Q 
. derived from QCD. 

^ ■ Our derivations in [l| included the one-loop mesonic corrections to the basic one-loop electromagnetic interactions, 
• w^ , so involved two loops overall. We derived a complete expression for the isospin breaking electromagnetic contributions 
^ ' to the baryon mass operator to this order, and then used the result to obtain expressions for the operators that lead 
M , to intramultiplet mass splittings. As indicated, our emphasis was on the loop expansion, rather than on the splittings 
- - 1 obtained in a perturbative expansion in powers of the chiral symmetry breaking parameters as in |9|| . The Feynman 
loop integrals that appear in the final results are of course the same as those obtained using the usual effective fields, 
but have been reduced in ^ for loops containing heavy baryons to integrals for time-ordered graphs. 

It remains to obtain explicit expressions for the mass splittings, calculate the integrals that appear, and compare 
the results with experiment, the subjects of the present paper. We will do this by using the natural connection of the 
"quark" description of the effective fields to semirelativistic dynamical models with the effective quarks playing the 
role of structure quarks. This will allow us to model the structure of the matrix elements beyond the point-baryon 
approximation and estimate the matrix elements of the effective Hamiltonian for finite-size baryons. The results 
will necessarily be model-dependent, but can be regarded as involving physically- motivated choices of cutoffs that 
emphasize the calculable long-distance parts of the loop corrections in the sense discussed by Donoghue, Holstein, 
and Borasoy JTd\. 

The electromagnetic contributions to the baryon masses can be expressed as a linear combination of the matrix 
elements of a set of independent spin- and fiavor-dependent operators F^, (i = 1,...,32) in the baryon states [30| . 
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There are 2 independent operators that involve only one set of quark indices, called one-body operators, 12 two-body 
operators, and 18 three-body operators. 

We showed in [l| that no three-body operators are generated by one-loop mesonic corrections to the initial two-body 
Coulomb and magnetic moment interactions after renormalization, even though diagrams that connect three sets of 
quark indices are initially present. The number of F's that can appear to this order is therefore reduced from 32 to 
14. When the calculation is restricted to the intramultiplet splittings, the number of independent matrix elements is 
further reduced to 4 [ll,[Il[ll. As a result, we can bring the electromagnetic mass-difference operator to the form 

niB^aJ2 + bTi + cFs + dTi3 , (1.1) 

i 

where 

the Qs are quark charge matrices, and the matrices Af^ and M" are defined as M"* — diag (0, 1, 0), M" = diag (0, 0, 1) 
[3l| . The complete set of one- and two- body electromagnetic operators is given in Appendix „ The procedures needed 
to convert Eq. ()l.ip to an expression in terms of effective baryon fields in the matrix representation are described in 
detail in [1, but will not be needed here. 

As will be discussed in Sec. HIl the results in [ij omitted two further isospin breaking contributions proportional to 
nid — rus that are allowed by the general effective field theory. These are analogous to the effective point interactions 
of the form ^ cr^ • (TjMf and ^ cr^ • a-jM^M^ that appeared in our earlier analysis of intermultiplet mass splittings 

0) Q , but involve the replacement of (one) matrix by M"* and have coefficients proportional to the mass difference 
(rud — m.u) rather than mg. As shown below, terms of this form arise in QCD-based dynamical models from a color- 
magnetic interaction [Ts^ , and can give IB contributions to the intramultiplet mass splittings comparable to those from 
electroma gnet ic interactions [l3, 111] [s^l • Since the extra effective interactions are two-body, the usual quark- model 
sum rules [1^ [TtI. [isj for baryon masses continue to hold for matrix elements. Once again, only four parameters a, b, 
c, and d are needed to describe the intramultiplet mass splittings through first order in (md — rriu), corresponding to 
an effective interaction of the form in Eq. (jl.ip . 

Our objectives in the present paper are, first, to extend the theoretical analysis of isospin breaking mass splittings 
by including the effects of the "color-magnetic" interactions [l^ missed in and, second, to estimate the final 
coefficients a, 6, c, and d in Tim numerically in the heavy baryon approximation. The first part of the analysis is 
completely general. The second will depend on specific models for baryon structure. 

As shown explicitly in Sec. Ill C[ the electromagnetic contributions to the coefficients a, 6, c, and d are given by 
combinations of cutoff-dependent Feynman integrals over the three momenta in pion and photon loops. The effect 
of the soft, long-range parts of the interactions are expected to be calculable. We will therefore make reasonable 
estimates of the soft parts of the integrals by employing known form factors and dynamical models in cutting off 
the divergent integrals. The differences from the experimental values give estimates of the sizes of the residual hard, 
short-range contributions. 

With the color-magnetic interaction included, we find that the calculated IB corrections are of the right general size 
and have the correct sign pattern to account for the pattern of the coefficients in the IB mass-difference Hamiltonian. 
The theoretical results agree with their experimental values at or within the margin of error, especially given the 
uncertainty in the theoretical input. Unless the corrections in the chiral expansion from two meson loops (three loops 
overall) are larger than we would expect, any deviations from the experiment can be attributed to the presence of 
short-distance effects which can be parametrized but not calculated in the effective field theory. 

II. ISOSPIN BREAKING CORRECTIONS TO BARYON MASSES 

It is necessary to recall that the leading IB corrections to baryon masses come from the Coulomb interaction, 
magnetic moments interactions, the effects of the d, u quark mass difference, and the color-magnetic interaction. We 
note that except for the color-magnetic interaction, we have included the one-loop mesonic corrections to other basic 
electromagnetic interactions (two loops overall) in our calculations. We wish to evaluate the contributions from these 
effects to the parameters a, b, c, and d. For this goal, we firstly discuss the contributions from the color-magnetic 
interaction to the mass splittings. 
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A. Color-magnetic contributions 

The QCD color-magnetic interaction was first introduced by Sakharov and Zel'dovich [l^ and further developed by 
De Riijula, Georgijand Glashow [l^ for their study on hadron masses in a gauge theory. Sakharov [iJl and Franklin 
and Lichtenberg [I^ showed that the QCD color-magnetic interaction was of the same order of magnitude as the 
purely electromagnetic interaction. The color-magnetic contributions to the mass of a baryon B arise in relativistic 
quark models from the interaction 
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(2.1) 



where rrii, rrij represent effective quark masses, and as is the strong coupling. At this point we will not attempt to 
calculate the matrix elements of this interaction, but will simply use its spin and mass structure to determine the 
equivalent effective operators to use with the heavy-baryon effective fields. We want, in particular, to see the extent 
to which the interaction in (Eq. (|2.ip ) can be reduced to an operator expression in terms of the F's and M''. For that 
purpose, we take to„ as the standard mass and consider the symmetrical case in which < d^{rij) >= A, a constant, 
for all states. We then substitute the identity [s^ 
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into Eq. (|2.ip and expand. To first order in — rriu, the color-magnetic interaction can then be rewritten as 

jB _ \ ' s/rd_ _ I \ ' n,fc!(i,fs_ _ , (2 3) 



where 



^ _ Anas A md-m^ ^, _ ^ m,, - m^, 



9 m?, 



(2.4) 



and the ellipsis represents terms that can be absorbed into the structures that already appear in our analysis of 
intermultiplet mass splittings in [5]. 

As mentioned above, the terms that are written out explicitly involve isospin-breaking operators of order — m„ 
that are allowed in the effective field theory [33| and could have been introduced from the beginning with unknown 
coefficients to be determined from the data. The derivation suggests that these represent short-distance effects not 
calculable in the chiral expansion. 

We now show how to put the color-magnetic interaction in the standard a, 6, c, d form. Using the identity 



we rewrite the operators in Eq. (|2.3p as 



-3g2 + -1 - 

3 



(2.5) 



^ Mt<T, ■(T, = -Y,cT,-cT,-Y, M/o", • a, - 6F2, 



(2.6) 



MfM'^a, • = ^ ^ Af/cr, • cr, - ^ Mf M;<t, • a, - GFn . (2.7) 

The first two operators on the right in each expression contribute to intermultiplet splittings but not to isospin 
splittings within multiplets, so these can be dropped. As far as the splittings are concerned, the following relations 
hold 

r2 = ^ E ^-^^ + ^(^4 - Ts) , Fn = -Fi3 + \{T^ ~ F5) . 
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Hence, we find 

"^CM — '^CM ^ + 6cMr4 + ccmTs + dcMri3 , (2.8) 
i 

where 

acM = -2C , 6cM = 3(C - C") , c = -3(C - C") , d = 6C" . (2.9) 
Note that the sum rules for baryon masses continue to hold [13, [lH ■ 

The one-loop mesonic corrections to the effective point interactions do not introduce new isospin breaking terms, 
are expected to be small, and will not be considered. 



B. Baryon mass splittings and sum rules 

The contributions of the operators Y^^ Mf, Fs, and Ti^ to the mass splittings within baryon multiplets can be 
determined from the results given in [U, |Tl[. Using those results (see, for example, Table I of the reference [l| [35j ) 
and Eq. (|l.ip . we find 
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(2.10) 



Since there are only four independent parameters in TYib, there are six sum rules among ten mass differences. Below 
are the well-known sum rules [ll [13, H, [Mi 

A" - A+ = n - p , 
A- - A++ = 3{n~p) , 

AO - A++ = 2{n -p) + (£0 - S+) - (S" - 2°) , 

S--SO = (S--S+)-(n-p), (2.11) 

S*--S*0-(I]*--S*+)-(n-p), 

2E*° - £*+ - S*" = 2i;0 - I]+ - S". 



These sum rules hold for any set of purely one- and two-body interactions as shown in [16|, |17|, UM, |21| . 

The sum rules can be violated by three-body operators. However, as shown in [1] and mentioned above, no effective 
three-body operators are generated through one loop in the chiral expansion, so any three-body effects must involve 
at least two meson loops and are expected to be small. The sum rules are therefore expected to hold with reasonable 
accuracy, as they do. 

A weighted fit to the seven known mass splittings other than those for the A baryons is given in Table H] A best 
fit is obtained at values (in MeV) of a = 1.82 ± 0.04, b = 3.35 ± 0.24, c = -1.78 ± 0.23, and d = 1.00 ± 1.40 with an 
average deviation from experiment of 0.13 MeV and a = 1-67 (with 7 degrees of freedom). Hereafter, we denote 
the electromagnetic mass-difference operator with the best-fit coefficients as TiJ'g'''. 

Using the data given in Table HI we find that there are no significant violations of the sum rules. 



5 



TABLE L A weighted fit to the seven accurately known baryon mass sphttings using the expressions in Eq. A best fit 

is obtained at values (in MeV) of a = 1.82 ±0.04, 6 = 3.35 ±0.24, c = -1.78 ±0.23, and d = 1.00 ± 1.40. The average deviation 
of the fit from experiment is 0.13 MeV. The experimental data are from [2^ . 



Splittings 


Calculated 


Experiment 


71 — p 


1.29 ± 0.12 


1.293 ± 0.000 


E- - E+ 


8.05 ± 0.62 


8.08 ± 0.08 


E- -E° 


4.81 ± 0.30 


4.807 ± 0.035 




6.75 ± 0.58 


6.48 ± 0.24 


E*" - E*+ 


4.48 ± 0.49 


4.40 ± 0.64 


E*- -E-» 


3.03 ± 0.33 


3.50 ± 1.12 


■=*- _ =*o 


3.19 ± 0.52 


3.20 ± 0.68 


A++ - A° 


-1.02 ± 0.56 




A++ - A- 


-3.88 ± 0.35 




A+ - AO 


-1.29 ± 0.13 





C. Expressions for the parameters a, b, c, and d 



We are now ready to determine the expressions for the parameters a, b, c, and d. Note that the IB mass-difference 
operator Hm can be written as jsH] 

HiB =Hi +H2 + ^-3 + ^CM ■ (2-12) 
The first term in this expression, Tii, is the total contribution to the baryon mass differences from charge interactions 



+24{i2,n - i2,K)ri3 - 2(Xi,^ - 2:i,K)(6ri3 - stio - 2ri4 + rn) 



(2.13) 



2X1, TT — 3^"^^'^ ~^2,k] 



[Mf + Mf + M^) , 



where 137 



12 2 2 



(2.14) 



In these expressions Iqq is an integral associated with the Coulomb interaction diagram in Fig. [U and li^i (i = 
1, ...,8; I — TT, K,r]) are the integrals associated with the diagrams shown in Figs. [H^a) -[Hc),[31IHa) -IH^d), respectively, 
that contribute to the baryon mass differences through charge interactions 38] . 

As discussed in [l[, we work in the heavy baryon limit, and the original Feynman integrals reduce to integrals over 
three momenta in old fashioned time-ordered perturbation theory. Thus, in Figs.[T][5l a solid vertical line represents a 
quark moving upwards toward later times, dashed lines represent mesons, wiggly lines represent transverse photons, 
and a horizontal dotted line represents the instantaneous Coulomb interaction between the particles on which it 
terminates. Only quark lines involved in the interactions and representative time orderings are shown. The incoming 
and outgoing quark lines are to be collected into the corresponding baryons so that, viewed at the baryon level of 
the usual baryon chiral perturbation theory. Figs. [1] and [3)Ja) are one- loop diagrams while Figs. [H [3Ib-d), and 2] are 
two-loop diagrams. Note that our Figs. [T][5] are identical to Figs. 2(b), 4, 5, 7, and 8 in [1], respectively. The evaluation 
of the various integrals is discussed in the next Section. 

The second term in Eq. (|2.12p , Ti.2 , is the total contribution to the baryon mass differences from magnetic moment 
interactions 



n2 



2 Ha fdtT 14) + n'g + n[^ 
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(2.15) 
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FIG. 1: One- loop electromagnetic corrections to the baryon mass due to the Coulomb interaction between quarks. Note that 
the "quarks" represent spin-flavor index sets on baryon flelds (not dynamical QCD quarks) and that spectator quarks are 
suppressed. 



FIG. 2: Two-loop corrections to electromagnetic interactions that involve meson exchange between quarks. 
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(c) 



Tla 



2 4 

-:^^^lI9,n{^2 + 3r4) + 6nl{l9,7r - -^19,7^)^5 

2 

+ 3Ma(-^9,7r - ^9,x)(9rio + Tn) + 4/1^ [/Ua/g.TT - (/"a + M6)-f9,if ] ^13 



(2.16) 
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FIG. 3: (a): The basic meson exchange diagram, (b), (c): 
electromagnetic correction to the meson-quark vertex. 



Electromagnetic contributions to the meson mass terms, (d); 




(a) 



(b) 



(c) 



FIG. 4; Mesonic corrections to electromagnetic vertices. 
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and 



(2.17) 



Here, = 2.793 , /it = —0.933; /^^ and /g ^ are the integrals associated with the direct interaction between magnetic 
moments (Fig. EJa)) and the moment-moment interaction including one-loop mesonic corrections (Figs. [5l^b) and 
[5]Jc)), respectively. 

The third term in Eq. (|2.12p . 7^3, involving the effects of the d, u quark mass differences on the baryon masses and 
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FIG. 5: Instantaneous magnetic moment-moment interactions and mesonic corrections. A zigzag line with crosses at the vertices 
represents a factor Ti^p = —fiifijCTi ■ (Tj /^^ 
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on the single meson exchange amphtude, is of the form 

Hs = ^'^t + ^Jxir^4,K" [i2ri3 ~ erio - ^v^ + 2rn] , (2.18) 

where Af = M|.o -M^± + A^4> ^em = ^'ll±-^'^lo, and Ad„ = [(md-mu)/(ms-m„)]a,„ is originally the coefficient 
of Mf from the single-particle mass operator at the quark level [s^. Using the value Um ~ 178 MeV obtained in 
the absence of loop corrections 0, [1] and the ratio (mj; — m„)/(ms — to„) k. 0.227 [1^ gives the estimate A^j^, w 4.04 
MeV. Here, we will take into account the possible loop corrections to a,„ and will treat Arf„ as a parameter. 

The operators Ti that appear above are all independent. However, contributions of their matrix elements to 
intramultiplet mass splittings satisfy a number of relations with Fio = — Tia, Fn = — — Fio -I- \iX4, — Ts), and 
T2 — — \ -^^i + 5(^4 ~ ^5)- Note also that Fig — F20 and F25 = F26 do not contribute. We can therefore bring 
Tij, {j ~ 1, 2, 3) to a form similar to Eq. (jl.ip 

= % + + CjFs + rfjTi3 , (2.19) 



where the coefficients of Hi are 



4 

bl = lQQ+4Ii^^+iIi,K-S{2l2,^+l2,K), 

ci = ~3Ii.K, (2.20) 
di ^ 12[-(Xi., -Xij<) + 2(X2,. -12,^)], 



the coefficients of H2 are 



2 

02 = gMa [(8/ia - Aifc)-?"9.7r + 2(^a + Mb)-f9,K] , 

^^2 = MaMb^W + ^Ma [(IS/la " l0fib)l9,7T " (IS^ia + 2'iHb)h,K " 12(/ia + Aifc)/9,,,] , 

C2 = -fJ.a{fJ.a + Mb)-^^ + ^A*a [10(Aia + Mb)^9,7r + (25/1^ + 2A^b)h,K + '^2{^la + Mb)^9,»;] , (2-21) 
4 

^2 = -2^Q^b/^^ + -/Xa [(5//;, - 14^a)^,7r + (8/la + 9Hb)l9,K + 6{fJ-a + /Ub)^9,J)] , 

and, finally, the coefficients of are 

a3 = Ad„, 63 = 2(Af /A,^4)/4,KO , C3 = -2(Af /A,^4)/4,KO , d3 = 8(Af /A,^4)/4,ko ■ (2.22) 
It follows from Eqs. (|2A2ll . and ((2l9| that 

3 3 3 3 

a = ^ -I- acM , 6 = ^ + bcu , c = ^ q -I- ccm , d = ^ + dcM ■ (2.23) 

i— 1 i—1 i—1 i—1 
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The coefficients acu, bcM, ccm, and dcM are defined by Eq. (|2.9p . 

In the next section, we will first evaluate the integrals and then study how well the dynamical theory developed in 
our earlier analyses of the baryon masses and magnetic moments describes the coefficients in Hib ■ 



III. THE INTEGRALS 



Let us start with the Coulomb integral Iqq and the magnetic integral /^^ defined as 

^^5" J (2^)3|fc|2' 3 J (2^)3' ^ ^ 

where /^^v = ^/{'^mN) is the nucleon magneton. 

Note that Eq. p.ip refers only to the "bare" Coulomb and magnetic interactions between quarks without effects 
of baryon structure. These integrals and others to follow are formally divergent. The divergences are commonly 
absorbed in the HBA to chiral effective field theory into unknown constants representing short-distance effects that 
are to be evaluated from experiment. 

In dynamical models, the integrals arise from matrix elements of the corresponding operators in physical baryon 
states, and those matrix elements involve additional momentum-dependent factors associated with the structure of the 
baryon which may regularize the integrals [40j . Our objective here is to estimate the soft or long-distance contributions 
to the integrals, and to see the extent to which these account for the magnitude and pattern of baryon mass splittings. 
We therefore adopt the dynamical approach. This requires information on the internal structure of the baryons, some 
of which is available through measured baryon form factors and the semirelativistic theory of baryon structure, but 
its use necessarily involves models that go beyond the HBA to chiral effective field theory. 

The semirelativistic theory of baryon structure has been considered by a number of authors and is quite successful 
[H, IHj S^- For simplicity, we will use the model considered in Q in which the baryon masses are calculated 
variationally for the semirelativistic Hamiltonian of Brambilla et al. using Gaussian wave functions. The results agree 
with those of a similar calculation by Carlson, Kogut, and Pandharipande [251 and are consistent with those of the 
much more extensive calculations of Capstick and Isgur [2^. 

We will use Jacobi coordinates to describe the positions of the quarks. Define 



■mi{xi - Xk) + mj{xj - Xk) 



ti-ijk - , (o-^) 

where the Xi are the particle coordinates, mij — nii + rrij, M = rrii + nij + mk, and Ryfc is the usual center-of-mass 
coordinate. The roles of i, j , and k are completely symmetric at this stage. However, it is reasonable to neglect the 
very small difference between the effective masses of the u and d quarks in the dynamical calculations. At least two 
of the quarks in each baryon are then identical or have the same mass. We label these 1 and 2, with the odd quark 
labelled 3 and then define the internal Jacobi coordinates p and A as p = ri2 and A = ri2,3. Alternatively, we can 
use coordinates with the role of (1, 2) replaced by (2, 3) or (3, 1) in the definition, and define p' = A' = r23,i, or 
p" = r3i. A" = r3i^2- The coordinate pairs p'. A' and p". A" can be expressed in terms of p and A and conversely, 
so one can work with whichever of the pairs is most convenient and switch between them as necessary. The spatial 
volume element is simply d'^Rd^pd^X, and is equivalent for the other pairs of internal coordinates. 
In the position space, we can express Iqq and 7^^, for the symmetrical case [ilj, as follows 

= ^V^'aM^, = ^/ "^'P^'^ mP,X)\'S^p) . (3.3) 

Using the position-space variational Gaussian wave functions in 7] for the L = ground states 

^o(p,A) = i^M^y\^p[-li(3y + (3lX^)] , (3.4) 
and changing coordinates appropriately 0, we can easily calculate Iqq and /^^. The results are 
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where acm = e^l^-K. For /3p = 340 MeV obtained for the nucleon, Iqq = 2.80 MeV and = 0.123 MeV. 
Next, we consider the integrals associated with the mesonic corrections. 

conies from the diagram in Fig. [2lja) and, as shown in [l|, factors into the product of a Coulomb integral and 
a mesonic integral ; . 



I : 



- 4/2 



•J2 



(3.6) 



where /3 is a common dynamical matrix element for meson emission Q [l^, Ei{k') = \/k"^ + M^-^, and Fyi(A; 2) is 
the axial vector form factor of the baryon introduced when we neglect the excited states and include the internal 
structure of the baryon through the baryon wave function. 

To evaluate the mesonic integral we use the dipole form factor FA{k — K\/{K\ + k'^Y used in our earlier 
analyses of baryon masses [1, [ij. The mesonic integral I[ ; is convergent and easy to be numerically calculated. For 
Ka = 850 MeV and (3 = 0.75, we find that = 0.039, I^k^ 0.016, and = 0.014. 

I2,i comes from the diagram in Fig. [2l^b) and is given by [43| 



hi = 



P 
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'2\\2 



(27r)3|fc| J {2Tr)32El E + El 



,{FQ{k')FA{k')) 



(3.7) 



Note that a product of the form factors, FQ{k^)FA{k is introduced at each vertex of the diagram to ensure that 
l2,i is convergent. To evaluate the integral numerically, we use a nucleon form factor FQ{k'^) — Aq/(Aq + fe^)^ with 
Aq = 843 MeV [13 . 

I^^i comes from the diagram in Fig. [H^c) whose intermediate matrix element, showing all quarks, is like a baryon- 
meson scattering matrix element. We will introduce a product of charge form factors FQ{k'^)FM{k'^) to the interme- 
diate state where FM{k'^) = A|j/(A^j + fc^)^ is the meson charge form factor with Am = 1017 MeV. Hence, 



l3,l = 



d^k 



4/2 / (27r)3|fc|2 



FQ{k^)FM{k^ 



l-^FA{k^)FA{k 2) 



{2^Y^E[E[' E[E, 



(3.8) 



where k" = k' + k . 

li^i comes from the diagram in Fig. [3lja) differentiated with respect to Mf . It is defined as 



4,i 



/32A 



M 



4/2 



/ {2^f2E[E[^-^^^^^^^ 



This integral can be analytically evaluated. The result is 

/32A,^4 \\{Mi + b\A) 
5127r/2 {Mi + XaT 



h,i = 



(3.9) 



(3.10) 



I^^i is the integral associated with the diagram in Fig. 4(a) and is identical to Ii^i, i.e., I5J = Ii^i. 

Iq^i comes from the diagram in Fig. [UJb). The extended structure at the vertex in Fig. [H^b) suggests that the same 
form factor as in Fig. [3]^ a) should be used. For the rest of the diagram, a Coulomb interaction must be absorbed by 
the wave function (not the form factor). So, for the symmetrical case, we get 



e2^2 
4/2 J (27r)3|fc|2 



d^k 



where 



^(fc,A) 



d''\ 



^(fc,A) 



d3fc' 



,Fi{k') 



(27r)32£;j 



d^p e 



-ik-p 



Using the Gaussian wave functions, Eq. ()3.4I) . we find 



^P{k,\) = exp[-A:V(4/3^)] , Iqq ^ e 



d^k 



(27r)3|fc|2 



(3.11) 



(3.12) 



(3.13) 
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TABLE II: Numerical values of U (i — 1, 9), Ii, and I2 for different meson loops. All the integrals are measured in MeV. 
Calculations use Qcm = 1/137, / ^ 93.0 MeV, (3 = 0.75, \a = 850 MeV, \q = 843 MeV, Xm = 1017 MeV, f3p = 340 MeV. 



Integral 


TT- meson 


K- meson 


7;-meson 


h 


0.108 


0.045 


0.040 


h 


0.113 


0.058 


0.053 


h 


0.059 


0.027 


0.024 




0.105 


0.024 


0.020 


h 


0.108 


0.045 


0.040 


h 


0.126 


0.087 


0.082 




0.100 


0.073 


0.069 


h 


0.057 


0.027 


0.024 


h 


0.005 


0.002 


0.002 




0.008 


0.025 


0.024 


I2 


0.038 


0.016 


0.014 



Hence, 



0^ 
4/2 



(2nY2E\ 



(3.14) 



Ij^i and Is.i come from the diagrams in Figs. ID^c) andlljd), respectively. The intermediate state of the diagram 
in Fig. md) has the baryon-meson scattering structure, thus it involves FQ{k'^)FM{k'^) ■ The diagram in Fig. is 
related to the one in Fig. IH^d) . Its intermediate part involves the baryon-meson scattering structure contracted with 
a meson-meson-baryon vertex but instantaneous, not scattering. Putting in the form factors, 



h.i - 



3^ 
Ap J (27r)3|fc| 



(fik' 



1 



-FA{k'^)FA{k"^)., 



and 



4/2 J (2^)3|fc| 



d^k' k' k'' 



{2n) 



E. 



-FA{k^)FA{k"'). 



(3.15) 



(3.16) 



At this point, it is important to note that the diagrams in Fig. [4] are related by electromagnetic current conservation. 
The sum of these diagrams associated with mesonic corrections to the photon-quark vertex must give a coefficient 
for the l/fc2 Coulomb singularity that vanishes in the limit of zero photon momentum, A; — )■ 0. This condition is 
found to hold for the diagrams without form factors [1]. It is straightforward to check that the condition still hold 
for the diagrams with wave functions and form factors. Indeed, since the integrals /g,; and I^^i appear with opposite- 
sign coefficients and for fc — > the integral over A in their expressions is just the normalization integral and hence 
approaches unity, the coefficient of l/fc^ vanishes when Iq i and are combined. Similarly, we can easily show the 
cancellation between the I^^i and terms for A; ^ if we write I^^i explicitly and notice that the form factors in 
the first factor in /g,; also approach unity for A; ^ 0. 

Finally, Igj and Iiqj come from Figs. [5jb) andl^Ic), respectively. They are identical and factor into the product of 
Ifj,^ and the mesonic integral I[ ; 



l9,i = I 



10,1 



— Ifi^i X I'l l 



(3.17) 



It is straightforward to evaluate Ii.i,l2,hl6,h and Ig^i numerically. To evaluate I^^i, Ijj, and Is,i, we integrate first 
on d^lk' with the polar axis chosen along k. The obtained results are then integrated numerically. 

We present in Table [TTl the numerical values of h {i — 1, 9), Ii, and I2 for different meson loops. All the integrals 
are measured in MeV. Calculations use aom = 1/137, / « 93.0 MeV, fi = 0.75, Xa = 850 MeV, Xq = 843 MeV, 
Xm = 1017 MeV, and /3p = 340 MeV. 

Since our loop calculations involve divergent integrals, cutoff-dependence of the results is inevitable. To explore 
this point further, we also estimate the integrals using a monopole form for the meson form factors suggested by the 
vector dominance model [2^. The cut-off parameters A's for the monopole form factors are chosen such that the mean 
square radii determined by the two forms (monopole and dipole) are the same. We show in Table IIIII the numerical 
values of /i (i = 1, 9), Xi, and T2 for different meson loops for case of using monopole form factors. All the integrals 
are measured in MeV. Calculations use = 601 MeV, Xq = 596 MeV, Xm = 719 MeV, acm = 1/137, / w 93.0 
MeV, /3 = 0.75, and /?„ = 340 MeV. 
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TABLE III: Numerical values of U {i — 1, ...,9), ^i, and X2 for different meson loops when using monopole form factors. All 
the integrals are measured in MeV. Calculations use Oom = 1/137, / f» 93.0 MeV, fJ = 0.75, fJp = 340 MeV, Xa = 601 MeV, 
Xq = 596 MeV, and Am = 719 MeV. 



Integral 


TT- meson 


K- meson 


77-meson 


h 


0.176 


0.097 


0.090 


h 


0.185 


0.110 


0.102 


To 


n 1 1 


n rifiQ 


U.UU J 


h 


0.132 


0.041 


0.035 


h 


0.176 


0.097 


0.090 


h 


0.194 


0.149 


0.144 




0.168 


0.135 


0.131 


h 


0.111 


0.069 


0.064 


I9 


0.008 


0.004 


0.004 


II 


0.052 


0.060 


0.058 


12 


0.045 


0.021 


0.020 



IV. THE COEFFICIENTS a, b, c, AND d 

To evaluate the coefficients a, 6, c, and d, we first need to estimate the coefficients acM, bcM, ccm, and dcM given 
by Eq. ((2^ . Using m„ fa 340 MeV, (to^ - m„) sa 2.5 MeV, {nid - m„)/TO„ = 0.0074 ± 0.0012, and a resuh from the 
simple quark model ^4J| 

6—^ A = mA - w 300 MeV , (4.1) 

we get C ~ 0.37 ± 0.06 MeV. Thus, for nis « 500 MeV, the color-magnetic contributions to the four parameters a, b, 
c, and d (in MeV) are 

acM = -0.74 ± 0.12, 6cm = 0.75 ± 0.12, ccm = -0.75 ± 0.12, dcu = 0.70 ± 0.44 . (4.2) 

Next, using = 1260 MeV^ A*^ = 5196 MeV^ and the values of the integrals given in Table Hi] (Table Hn] for 
the case of using monopole form factors), it is straightforward to calculate numerically the coefficients b, c, and d. To 
evaluate the coefficient o, we have used A^u as a parameter and fitted a exactly to its experimental value. A best fit 
is obtained at Adu — 2.54 ± 0.00 MeV (Adu — 2.43 ± 0.00 MeV when using mono pol e form factors), a value that is 
lower than the value A^^ « 4.04 MeV estimated ignoring loop corrections to ^M- The calculated values of the 
coefficients a, b, c, and d (in MeV) are 

a = 1.82 ± 0.12, b = 3.00 ± 0.12, c = -1.46 ± 0.12, d = 2.28 ± 0.44 , (4.3) 

for the case of using dipole form factors, and 

a = 1.82 ± 0.12, b = 3.34 ± 0.12, c = -1.50 ± 0.12, d = 2.48 ± 0.44 , (4.4) 

for the case of using monopole form factors. 

The calculated values of a, 6, c, and d for Hm and the contributions to those values from Hi, H2, W3, and Hcm 
for these two cases are shown in Table HVl The best fit values of the parameters from Tijg*'* are also given there. 

The results given in Table IIVI show that the calculated IB corrections are of the right general size and have the 
correct sign pattern. These corrections seem to account fairly well for the pattern of coefficients a, 6, c in the IB mass- 
difference Hamiltonian. The coefficient d has the correct sign, but its calculated magnitude appears to be somewhat 
large. However, like the other calculated values of a, b, and c, the calculated value of d agrees with its experimental 
value within the margin of error. Using the values of a, 6, c, and d given in Eq. (j4.3p for the dipole case, we can easily 
evaluate the baryon mass splittings shown in Eq. (|2.10p and find that the average deviation of the calculated values 
from experiment is 0.23 MeV with a weighted /n/ of 1.63 for n/ = 3 degrees of freedom {p « 0.2). 

Note that employing a monopole form for the meson form factors improve the results for the coefficients b and c, 
but not for d. Nevertheless, the calculated value of d still agrees with its experimental value within the margin of 
error. Using the values a, b, c, and d given in Eq. (|4.4p for the monopole case, we find the average deviation of 
the calculated values from experiment to be 0.29 MeV with a weighted X^/'^/ = 0-86 for nj — 3 degrees of freedom 
(p « 0.36). The weighted fit is clearly better even though the average deviation is larger, the result of the uncertainty 
in the experimental value of d. 
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TABLE IV: Coefficients a, b, c, and d from 7ii, TI2, Ti-s, TCcm, Ti-m, and Ti^g'^'. All values of the coefficients are in MeV. 
Calculations use A^^, = 1260 MeV'^, A^^ = 5196 MeV^, the calculated values of the color-magnetic parameters given in Eq. 
(1421) . Adu = 2.54 ± 0.00 MeV obtained from a fit, and the values of the integrals shown in Table |ll] (Ad„ = 2.43 ± 0.00 MeV 
and Table Hm for the case of using monopole form factors). The uncertainties shown for the theoretical values of the coefficients 
include only those from the uncertainties in the color-magnetic terms. 



Hamiltonian 




Dipole 


case 




Monopole case 




a 


b 


c 


d 


a 


b 


c 


d 


Hi 


-0.05 


2.17 


-0.08 


0.74 


0.01 


2.29 


-0.18 


0.67 


H2 


0.07 


-0.12 


-0.43 


0.04 


0.12 


-0.03 


-0.23 


-0.23 


Ti.3 


2.54 


0.20 


-0.20 


0.80 


2.43 


0.33 


-0.33 


1.34 


'Hem 


-0.74 


0.75 


-0.75 


0.70 


-0.74 


0.75 


-0.75 


0.70 


HlB 


1.82 ± 0.12 


3.00 ± 0.12 - 


1.46 ± 0.12 


2.28 ± 0.44 


1.82 ± 0.12 


3.34 ± 0.12 


-1.50 ± 0.12 


2.48 ± 0.44 


T/bcst 


1.82 ± 0.04 


3.35 ± 0.24 - 


1.78 ± 0.23 


1.00 ± 1.40 


1.82 ± 0.04 


3.35 ± 0.24 


-1.78 ± 0.23 


1.00 ± 1.40 



Since most of the integrals clearly have potentially large short-distance contributions (the divergences in the absence 
of form factors are strong) , the presence of the short-distance effects is likely the source of any deviations from the 
experiment results. We do not regard the distinction between the dipole and monopole fits as definitive given the 
different short-distance behavior of the corresponding integrals. 

To illustrate the uncertainty within the model, we note that if one supposes that the quark wave functions are 
exponential at short distances rather than Gaussian, exp (— r/a) rather than exp {—0j^r'^/2), but keep the dependence 

on A and the mean square separation < > the same, then a — l/(-\/2/3p) and the Coulomb integrals increase by 
^-K 12. That is, the Coulomb integrals are underestimated using Gaussian wave functions which neglect short distance 
correlations. The more complicated wave functions used by Carlson tt al. [25^ include correlations and indeed give a 
some what larger energy as noted in ^Ij. The same remarks hold for the magnetic integral. The short distance effects 
are even stronger there: the magnetic integrals calculated using the exponential wave functions increase by a factor 
2\ph^. Therefore, Iqq and /^^ can be treated as parameters by multiplying their expressions defined in Eq. (j3.5p 
with the scale factors. We find, however, that doing so does not substantially improve the results. 

It is encouraging that the model, with all parameters except A^u fixed gives as an accurate account of the purely 
electromagnetic parameters 5, c, d. 

V. CONCLUSION 

Our results here consist of a numerical analysis of the IB contributions to the baryon masses which we analyzed 
recently using a loop expansion in the heavy baryon approximation to chiral effective field theory and methods 
developed in our earlier analyses of the baryon masses and magnetic moments. The leading IB corrections to the 
baryon masses come from the Coulomb interaction, magnetic moments interactions, the effects of the d, u quark mass 
difference, and effective point interactions attributable to color-magnetic effects. 

We have made reasonable estimates of the various integrals by introducing the known form factors and wave 
functions from successful semirelativistic models for the baryons and using the results to evaluate the parameters. We 
find that the resulting contributions to the IB corrections arc of the right general size, have the correct sign pattern, 
and agree with the experimental values within the margin of error. We also find that employing monopole instead 
of dipole form factors slightly improve the results, but scaling the Coulomb or magnetic moment contributions, does 
not. To the extent to which effects of adding a second meson loop are small, a view supported by the smallness of the 
three-body violations of the sum rules, it appears likely that any deviations from the experiment can be attributed 
to the presence of short-distance effects which can be parametrized but not calculated in the effective field theory. 

APPENDIX: ONE- AND TWO-BODY OPERATORS 

We present here sets of one- and two- body operators defined earlier in [H, Q. In the case of O(e^) contributions to 
the baryon masses, the F's must be bilinear in the quark charge matrix Q = diag (2/3, —1/3, —1/3) and can depend 
otherwise on the quark spin matrices cr and and fiavors. Ignoring isospin breaking through the small u, d mass 
difference and using the conventions that 

Em-E' E[^^i-^E. Efo-^]-^ E' 
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where k G u,d, s label the three quarks in a baryon, we can group the F's into sets of one- and two-quark operators 
as follows: 

One-body operators: 

ri = ^[Q?], ry^^[Q,?M/]. (A.2) 

Two-body operators: 

T2 = J2iQU<^.-(^j)i r4 = ^[Q,Q,], r, = Y.^Q,Qj{(T,-rT^)], (A.3) 

Ts = J2[QlM:ia,-a,)], Tio-^iQ^M;], Tn - ^[Q?A/;(^. • <t,)], (A.4) 

ri3 = ^[Q.Q.M/], = ^[Q,Q,M/(rT, • ^,)], (A.5) 

ri9 = 5][Q?MfM;], = ^[q?m/m;(<t, • (a.6) 
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